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foreground structure? 
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ABSTRACT 

The steep source counts and negative iGcorrections of bright submillimetre galaxies 
(SMGs) suggest that a significant fraction of those observed at high flux densities may 
be gravitationally lensed, and that the lensing objects may often lie at redshifts above 
1, where clusters of galaxies are difficult to detect through other means. In this case 
follow-up of bright SMGs may be used to identify dense structures along the line of 
sight. Here we investigate the probability for SMGs to experience strong lensing, using 
the latest iV-body simulations and observed source flux and redshift distributions. We 
find that almost all high redshift sources with a flux density above 100 mJy will be 
lensed, if they are not relatively local galaxies. We also give estimates of the fraction 
of sources experiencing strong lensing as a function of observed flux density. This 
has implications for planning follow-up observations for bright SMGs discovered in 
future surveys with SCUBA-2 and other instruments. The largest uncertainty in these 
calculations is the maximum allowed lensing amplification, which is dominated by the 
presently unknown spatial extent of SMGs. 

Key words: submillimetre - gravitational lensing - galaxies: statistics - galaxies: 
high-redshift - cosmology: observations. 



1 INTRODUCTION 

Sources that are picked up in deep extragalactic millimetre 
(mm) and submillimetre (submm) wavelength surveys are 
far from being ordinary galaxies. Due to strong evolution, 
and also driven by the typically coarse angular resolution of 
current surveys at these wavelengths, mm/submm sources 
are rare and luminous compared with well-studied popula- 
tions of optical galaxies. Submm galaxies (SMGs) are usu- 
ally interpreted as an early rapidly star-forming phase, per- 
haps driven by a major merger, in the sequence t hat will ulti- 
mate ly become a massive elliptical galaxy (e.g.. [Blain et all 
|2002| . and references therein). 

However, gravitational lensing can also play a role, 
such that some fraction of SMGs may be more typical, 
lower-luminosity galaxies at high redshift (z > 1) that are 
strongly amplified by lensing. This is already a well known 
effect for quasars, where the steepening source counts make 
lensin g of increasing importance for the brightest objects 
(e.g., lKochanekll2006l ). In addition, the prevalence of gi- 
ant optical arc s is partly explaine d by the high redshifts 
of the sources i|Wu fc Chiuehfl2006D . There are many anec- 
dotal examples of SMGs that are either known or strongly 
suspected to be lensed. Examples include sources seen 
through targete d rich galaxy clus t ers, most spectacularl y 
in Abell 2218 (|Kneib et all Il996l ; ISwinbank et~ai1 120031 ). 



and MS0451-03 (jBorvs et alJl2004h . as well as the bright- 
est source in the central Hubble Deep Field North, usually 
called HDF850.1 (|Dunlop et alj|2004 ). the brightest source 
in the entire GOODS-North field, referred to as GN20 (Pope 
et al., submitted), and oth er bright well-studied sources, 
such as SMMJ14011+0252 dlvison et alj|200l l: [Smail et all 
120051 ). MIPS J142 824.0+352619 JBorvs et ail I2OO6I ). and 
SMM J02399- 0136 dlvison et al.lll998l ). It has long been rec- 
ognized (e.g., iBlain fc Longaii Il993l ) that the steep source 
counts and high redshifts of SMGs make the brightest ones 
prime lensed candidates. 

Despite some efforts to m odel the lensing contribu- 
tion to the source c ounts (e.g., IPerrotta et al ] 12001 120031 : 
iNegrello et alj |2007j ) the predictions for SMG lensing are 
still very uncertain. The purpose of the present paper is to 
improve this situation, using the latest numerical models, 
and to try to quantify the uncertainties. In the process of 
this investigation we will also be able to address a related 
question - namely whether bright SMGs could be used as 
tracers of galaxy clusters at redshifts that are high enough 
to be challenging for other techniques. 

In order to carry out this study we will need several in- 
gredients, including the unlensed submm source counts and 
the redshift distribution. However, we start in section [2] by 
considering the probability distribution for lensing ampli- 
fication along random lines of sight. We then look at the 
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observed and modelled source counts for SMGs in section [3] 
and use these to calculate the source counts after lensing in 
section [4] also exploring the dependence on the redshift dis- 
tribution of sources. A major simplification of our approach 
is to assume that the shape of the source counts is indepen- 
dent of redshift; in order to test this we calculate in section^ 
lensing expectations from a specific evolutionary model that 
is consistent with a range of current IR and sub-mm data. 
Finally, we consider the uncertainties in the various compo- 
nents of our analysis in section [5] and interpret our results 
in the context of upcoming submm surveys in section [7] 



2 PROBABILITY DENSITY FUNCTION 
2.1 Millennium Simulation 

The first ingredient we will need for our stud y is the sta- 
tistica l distribution of lensing amplifications. IWang et al.l 
(12002) 1 presented a 'universal probability distribution func- 
tion' valid for all cosmologies and redshifts. However, it is 
only applicable to weak lensing, with a maximum magni- 
fication much less than 10. As we are interested in strong 
lensing effects leading to the very brightest sources, where 
we will find a mean magnification far beyond the weak lens- 
ing regime, t his is unsuitable for our purposes. Other rele- 
vant studies (|Perrotta et alj|2002i ; iNegrello et al.ll2007f l have 
used largely a nalytic approximations, have also fo cused on 
weak lensing IjJain et al.l 120001 ; iKeeton et al. 200 5]), or have 
relied on TV-body simulations (|Bartelmann fc W eiss 1994; 
IWambsgansl Il992) which have since been dwarfed in size 
by the Millennium Simulation l|Springel et al.ll2005l ), a com- 
putational run containing approximately 10 10 particles of 
mas s 8.6 x 10 s Ma/fc in a cube of side 500 Mpc/ft. 

iHilbert et al.I (|2007) studied lensing amplification prob- 
ability distributions using a ray tracing approach with the 
Millennium Simulation, and a ACDM cosmology with pa- 
rameters Qm = 0.25, Qa = 0.75, h = 0.73, n = 1, and 
as — 0.9. They calculate optical depths and for an im- 
age having been magnified by a factor fi. The former gives 
the optical depth assuming a uniform distribution of images, 
whereas the latter assumes a uniform distribution of sources. 
Since we are lensing an approximately random background 
of submm galaxies, henceforth we will use only and its 
related probability density function (PDF) p s (/i), omitting 
the superscript 'S' from now on. 

Notice that this is not strictly the probability for a 
source undergoing a total magnification fj,, but only for 
one image of a source. There is some ambiguity here over 
whether we should be considering the total amplification for 
all images - this will depend on how SMGs are detected, the 
resolution of the imaging and on details of individual lensing 
models. Having pointed out this ambiguity, we will simply 
ignore it for the rest of the paper. 

Given p(fi) , the probability that a source has been am- 
plified by a factor /i in the range (fi, /i + d/x) is then p(fj,)d[i. 
Though the scalar magnification /i has a sign associated with 
it, where negative valu es indicate a change in image polarity 
(for a review see, e.g., ISchneider et al. 1 1 19921 ; ICourbin et al.l 
2002), the behaviour of the PDF at fi < 1 will allow us to 
take /i > 0. We will require that the PDF be normalised: 



p(fi)dn = 1. 



(1) 



Since flux of a source must be conserved, sources amplified 
in some directions must be deamplified in others. This flux 
conservation constraint takes the form (Hil bert et alj 2007) 



\n\p(n)dn = 1. 



(2) 



The ray tracing techniques used essentially assume 
point sources, which allows for arbitrarily high magnifi- 
cations. The consideration of extended obj ects imposes a 
maxim um magnification limit, estimated by IPerrotta et all 
(2002) to be in the range 10-30 for 1-10 kpc sources at red- 
shifts 1-4. The applicability of this estimate to SMGs, how- 
ever, is uncertain, since the characteristic size of emission is 
unknown. SMGs appear t o be extended objects of a few kpc, 
but m ay well be clumpy (|Pope et al.ll2005l ; IChapman et alj 
2003), so small knots of bright star-formation would have a 
higher potential for amplification by lensing than the galaxy 
considered as a whole. Perhaps th e best empirical constraint 
comes from I Younger et al.l (|2007l ). who recently constrained 
the sizes of the submm emitting regions to < 10 kpc using 
interferometry. Until we have much more concrete informa- 
tion about the size distribution of SMGs there is little we 
can do about the effects of finite source size. For now we 
continue without a maximum magnification limit, but will 
consider in section lBTTI how imposing such a constraint would 
affect our results. 



2.2 Transform to other redshifts 



IHilbert et all i|2007h give numerical data for a source plane 
at z = 2.1, as well as the peak position At p (z) and full width 
at half maximum FWHM(z) of the distributions for other 
redshifts. We found that we could transform the z = 2.1 
PDF to other redshifts using a transformation of the form 

V > A 4 = I ~ZhT7, H^p + ^Mp, (3) 



p(m) 



I gi>/M _ 1 

Ap(ji), 



(4) 



where /i p is the peak amplification at the reference z, 8/j.p 
is the shift in the peak between redshifts, b is determined 
by requiring that we obtain the correct FWHM, and the 
normalisation constant A is determined by requiring that 
we satisfy equation ((T}. The peak shifts to lower as z 
increases, to balance the high /i tail, allowing the PDF to 
satisfy equation ([2]). 

This transformation preserves the asymptotic be- 
haviour of the PDF, as discussed in section 12.31 which is 
most crucial to our analysis. We plot the resulting family of 
curves in Fig. \T\ For reference the values of the coefficients 
are (A, b) = (0.473,-1.504), (0.968,-0.051), (1.383,0.655), 
(1.681,1.077) and (1.889,1.346) for z = 1, 2, 3, 4 and 5, 
respectively. 



2.3 Limiting behaviour 

The PDF f or z = 2.1 is seen to h ave a fi~ 3 tail, as predicted 
by theory (Schneid er et al.lll992l ) above an amplification of 
about /ihigh — 20. Thus, the data are fit to 
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from the SCUBA Half Degree Survey (SHADES) at 850 fim. 
They use both a broken power-law of the form 
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Figure 1. Family of lensing amplifica tion PDF cur ves at different 
source redshifts, derived from Hilbcrt ct al. (2007). 



p{p) = ^highA* for fi > ^high, 



(5) 



where Ahigh is an adjustable parameter, and this functional 
form can be used to extend the PDF to arbitrarily high am- 
plifications. The /Xhigh limit is adjusted for all other redshifts 
according to the transformation described above. Inaccuracy 
in the fi~ 3 behaviour above this limit may cause a small kink 
in the function at the transition between the data and the 
extrapolati on, but we find this to be quite small. 

In the iHilbert et ail (|2007h data for z = 2.1 there is a 
cut-off at about fi = 0.8, where the PDF stops decreasing 
and reaches a near constant plateau. In this paper we are 
not interested in the deamplification effects from the fi < 1 
regime, so at most this low-/x behaviour is important only for 
the normalisation of the PDF. We find that by instead allow- 
ing the PDF to drop off like a Gaussian the normalisation 
of the PDF is changed by about 10 -6 , and so is negligible 
for our purposes. We therefore assume that the PDF goes 
to zero like a Gaussian for small fi, i.e., 



p(fi) = A lov , exp - 



- Mo) 2 



2<r 2 



for pb < m 

o 



(6) 



where Ai ow , /io, and a are adjustable parameters, and here 
ziiow is chosen to make as smooth a transition from the data 
to the Gaussian as possible. The peak of the distribution 
itself remains defined by the data. For example, for z = 2.1, 
the peak of the distribution is at /i ~ 0.95, while the low 

Gaussian tail is not used until /i < Mi ow = 0-8. 

This approach is in contrast to IPerrotta et alj (|2003l ) 
who approximated the entire weak lensing regime below /i ~ 
1.5-2.0 as a Gaussian and the entire strong lensing regime 
above the same limit by the pT A power-law. Our approach is 
to use numerical data for a broad intermediate range, giving 
a more accurate shape for the PDF, particularly for small 
and intermediate amplifications. 



3 SOURCE COUNTS 

3.1 Fits to observed counts 

Functional fits to observed diff erential source coun ts of sub- 
millimetre galaxies are given in lCoppin et al.l (|2006l ). derived 



diV 
~dS 



N' 
S> 
N' 
S> 



if S < S', 

if s > s', 

and a Schechter functional form 
diV N' ( S \ a+1 _ s /s> 



?' \S'J 



(7) 



(8) 



dS S> 

with units of mjy _1 deg~ 2 . The differential counts dN/dS 
can be integrated from S to infinity to obtain the cumulative 
counts N(> S) with units of deg~ . Currently there are not 
enough high flux (S > 20 mjy) sources to prefer one function 
over the other. 

We have chosen to use the Schechter function as our 
primary estimate of the source counts. The power-law drops 
off much slower than the exponential part of the Schechter 
function, allowing for many more intrinsically bright objects. 
As such, the most gradual power-law consistent with the 
data can be used as an upper estimate to the counts, which 
results in a far more conservative estimate of the strong 
lensing probabilities. A lower limit to the source counts is 
less helpful, as the corresponding uncertainty in the upper 
limit of lensing probabilities is unknown. 

The Sch echter function param eters that we adopt are, in 
parallel with lCoppin et all (j2006l ). N' = 1600, 5" = 3.3, and 
a — —2. For the maximal broken power-law we use N' = 58, 
S" = 9, q = -2.2, and (3 = -4.2. The contribution to the 
overall background (i.e. J S(dN/dS)dS) from a power-law 
number counts distribution diverges at the faint end if the 
slope a ^ — 2 (in equation (J7J)). To prevent overproducing the 
submm background, we therefore set a limit on the broken 
power-law of S ^ 0.2 mjy. 

3.2 Phenomenological models 

lLagache et ail (|2003h used a phenomenological approach to 
model galaxy evolution and predict number counts at a va- 
riety of wavelengths, including 850 /im, an d corr ected this to 
include constraints from Spitzer data in |2004|. Thei r stan - 
dard model is roughly consistent with ICoppin et alj (l2006l ) 
in the range where data from SHADES exist, but decreases 
much more slowly above about 50 mjy. This is due entirely 
to local bright Euclide an counts, for which dN/dS oc 5 ,_5/ ' 2 . 
I Chary fc E lbaz (2001) similarly have a model for the source 
counts that has been adjusted to include the Spitzer data. 
The results o f their model at 850 fim is very similar to 
lLagache et al.l (2004), including the bright Euclidean part 
as well. 

Normalising a Euclidean source count estimate to 
the SCUBA Loca l Universe Galaxy Survey (SLUGS; 
iDunne et alj |2000| ) . we find agreement wi t h the high 
flux behaviour in b oth the lLagache et al.l |2004) and 
ICharv fc Elbaj l|200lj ) models. Since the Euclidean part of 
the source counts comes from low redshift sources, we do not 
expect there to be any significant lensing here. Subtracting 
the S~ 5 ^ 2 part from either of these models to estimate the 
high redshift population, we find that the high flux fall-off is 
very similar to a Schechter function. This is only a first-order 
approximation, as in principle the source counts change for 
each redshift slice. We will investigate redshift dependence 
in section |4~21 and the change in the shape of the counts with 
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Figure 2. Differential number counts prior to applying tens- 
ing. S olid and long-dashed lines are the models of lLagache et all 
(|2004l l and lCharv k Elbaj [120011) . respectively. Short-dashed and 
dotted lin es are the Schechter functional form and broken power- 
laws from lCoppin et al.l ||2006|) . 



redshift in 15.21 We will also compare with the results from 
an evolutionary model in section [U But for now we will as- 
sume a single screen of sources with counts independent of 
redshift. 

Between the low flux limit of the models at about 1 /zJy 
and the knee at several mjy, both models predict inter- 
mediate counts between the Schechter and broken power- 
law fits. Fig. [2] compares these functional fits to the mod- 
els of lLagache et all (|2004l ) and lCharv fc Elbad (|200lh . This 
justifies our use of the Schechter and power-law functional 
forms as our 'best' and 'maximal' estimates, respectively. 
The lensed source count estimates in this paper at high flux 
densities should be considered to be in excess of known lo- 
cal (e.g., IRAS) galaxies, unless otherwise stated. We will 
see that, as expected, highly amplified sources make an in- 
creasing contribution at brighter fluxes. 



4 CALCULATING LENSING EFFECTS 
4.1 Observed source counts 

If an object with flux density S is amplified through lens- 
ing by a fa ctor /j, the object is observed with flux density 
Sobs = nS (|Schneider et alj|l992h . For a given observed flux 
Sobs, we run through a range of magnifications, find the cor- 
responding S, and assign p(fi)dN/dS to be the contribution 
to the lensed differential source counts from that /j,. In this 
way we can express the differential source counts after lens- 
ing as 

dN [°° AN ( S ohs \ , , , 

For different values of Sobs we plot the contribution to 
the above integral as a function of \i in Fig. [3] We do this 
specifically for sources at * = 2, as well as showing a separate 
panel with results normalised by the total contribution from 
all (i. In the normalised case, the line can be interpreted as 
a local p(jtt) for that Sobs- Although we have only plotted 
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Figure 3. Magnification distributions. The contribution to 
dAydSobs made by each fi for various Sobs at 2 = 2 (top) and 
normalised by the total dN/dS b s (bottom). The bottom panel 
shows how the high fi tail dominates over the fj, ~ 1 peak for 
all S > 30 mjy. The additional thick line in the top panel is the 
'maximal' estimate for Sobs = 10 mjy (the difference is typical 
for all Sobs)- 

results for z = 2, other redshift choices look qualitatively 
similar. 

We can also examine the contribution to the source 
counts coming from various magnifications by specifying a 
minimum \x that we might be interested in. The fractional 
contribution to diV/dSobs from /j, > pi m i n is found simply by 
changing the lower limit on the integral in equation [9] from 
to /imin, and normalising by the full dN/dS b s - This is 
plotted in Fig. 3] It can be seen that at high flux densities a 
large portion of sources will have significant amplifications. 

4.2 Comparing Redshifts 

Lensing PDFs and a dN/dS function can be combined at 
arbitrary redshifts in the range 0.5 < 2 < 5.7 (where 
iHilbert et alj (|2007l ) quantify their lensing simulation). The 
differential source counts of the lensed populations are found 
via equation [5] and shown in Fig. [S] for a range of source 
redshifts. One can see that the lensing tail is relatively im- 
portant above about 20 mjy and more so for higher redshift 
submm galaxies. 

The cumulative source counts, 

N(>S)= / — — dSobs (10) 

J s QOobs 

are used to determine the density of sources above a given 
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Figure 4. Combined contribution to djV/dS b s from all magni- 
fications above £t m i n for sources at z = 2. 
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Figure 5. Lensed differential source counts found using equation 
[9] Dashed lines represent counts derived from the best Schcchtcr 
function (lower solid line) and dotted lines those from the maxi- 
mal power-law (upper solid line). Each set of curves is for sources 
placed at z = 1, 2, 3, 4, and 5 from bottom to top. 



flux on the sky. This is a useful quantity for making obser- 
vational predictions for surveys. 

Examining the contribution to source counts as a func- 
tion of fimin, as described abov e, we set, somewhat ar- 
bitrarily (though co nsistent with iNegrello et al. 1 120071 and 
IPerrotta et"all l2002 ) , a characteristic magnification for lens- 
ing events of /i m i n = 2. In Fig.[6]we plot the contribution as a 
function of S^bs for various redshifts, as well as the 'best' and 
'maximal' estimates for the average over all redshifts. This 
quantifies the probability that an object observed at some 
S bs has been lensed with a magnification greater than 2. 
One could easily obtain similar results for any other choice 
of minimal amplification. 

There have been several attempts to characterise the 
redshift distribution of SM Gs. Though phenomenological 
models do exist (e.g.. iGranato et al.ll200 J ) we focus on esti- 
mates based purely on observational data. IChapman et al.l 
(2005) obtained a large number of spectroscop i c red shifts 
for radio selected galaxies, while iPope et alT (|2006h ob- 
tained both spectroscopic and photometric redshifts in the 
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Figure 6. The fraction of dN/dS oos due to magnifications greater 
than 2, which can be interpreted as the probability of a lensing 
event as a function of S ODS - The solid lines use the Schechter 
source counts and dashed lines use the maximal broken power- 
law, resulting in a minimal estimate of the lensing probability. 
We show (from bottom to top) z = 1, 2, 3, 4, and 5 as thin lines 
and the average over all redshifts as a thick line in each case. 



GOODS-North survey, and lAretxaga et all l|2007h used ra- 
dio/submm photometric redshifts in the SCUBA Half De- 
gree Survey (SHADES). All three are consistent with a 
Gaussian distribution centred at z — 2.2 ±0.1 and with 
a = 0.8 ± 0.1. This is the distribution we use to obtain the 
average over redshift in Fig. [6] The results are nearly iden- 
tical when varying the peak redshift by ±0.1. Fig. HJ] shows 
that for sources with S ODS — 20 mjy the probability of strong 
lensing is already 10-20 per cent for any reasonable redshift 
distribution. Empirical evidence suggests that this may be 
a conservative estimate, since at the moment all sources de- 
tected with Sobs > 20 mjy at 850 /im have been claimed to 
be lensed. 

We can also look at the average amplification as a func- 
tion of observed flux density, given by 



<M> = 



(ii) 



This is shown in Fig. [7] across a range of redshifts. One can 
see a change in behaviour at higher and lower flux densi- 
ties, depending on the source redshift. Below about 35mJy 
the mean amplification increases with redshift, while above 
35mJy it decreases with redshift. This is due to the change 
in dominance between the fi ~ 1 peak and the high \i tail, 
an effect that can be seen in the bottom panel of Fig. [3] 



5 UNCERTAINTIES 
5.1 Lensing probabilities 

iHilbert et af] ||2007h discuss how they may have underesti- 
mated the lensing properties by about 15 per cent, based on 
replacing the principle len sing objects in their ray t racing 
by analytical haloes of the iNavarro, Frenk fc White] (|l997T ) 
form. This is surely a conservative estimate of the effects 
since baryons are neglected. Most large-scale simulations, 
like the Millennium Simulation, include only collisionless 
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Figure 7. Mean magnification as a function of S ^ a . 
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Figure 8. Differential source counts including the effects of lcns- 
ing, but with a hard upper limit on the amplification of fi < 30. 
Here we have plotted lensed counts for the Schechter (dashed) 
and power-law (dotted) models, for the average over all source 
redshifts. This should be compared with Fig. \E\ making it clear 
how dramatic the effects are at the bright end. 



dark matter and not the much more complicated dissipa- 
tive baryons. Clum ping due to baryons is b ound to increase 
strong lensing fe.g.. lBartelmann et al.|[l995l ) and so the nor- 
malisation of the high /i tails that we us e should be consid- 
ered a s a lower limit. Recen t studies fe.g. JWambsganss et al.l 
|200S| ; iHilbert et all 120081 ) suggest that structure due to 
baryon cooling and dissipation increase the incidence of 
strong lensing by at least 25 per cent. 

However, the effect of finite source size goes in the oppo- 
site direction. For extended sources there is an upper limit 
on the magnification, as mentioned in section[2] The effect of 
this on the lensing probabilities can be significant, depend- 
ing on the counts and on details of the size distribution. We 
can use Fig. [7] to gauge how important this might be. If we 
expect the limit to be ~ 30, then the effect will be negligible 
for sources with Sobs J$ 30mJy, but increasingly important 
for brighter sources. 

We further show the importance of such a cut-off in 
Fig. [9] where we have calculated the differential source 
counts using the lensing PDF, but not allowing any events 
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Figure 9. Fraction of the differential counts having amplification 
H > 2, but restricted to lensing distributions with maximum of 
li = 30. Here we plot the results for the Schechter (solid) and 
power-law (dashed) counts models for the average over all source 
redshifts. This plot should be compared with Fig. [6] 



at /i > 30. One can see significant changes for the brightest 
counts compared with Fig. [5] We also show the changes in 
the probability of strong lensing as a function of observed 
flux density in Fig. [6] again using a maximal amplification 
of /i = 30. One can see that the fraction of strongly lensed 
sources is significantly reduced, particular for the power-law 
counts model. 

Is 30 a reasonable limit for lensing amplification? One 
indication that there is not a hard cut-off in fi is that the 
SMG lensed by the cluster Abell 2218 ha s been estimated t o 
have a (combined) amplification of ~45 |Kneib et al.ll20o3 ). 
However, since there is currently only vague size informa- 
tion for the submm-bright regions of SMGs, there are no 
strong constraints on the maximum amplifications that are 
possible, nor on the shape of any distribution at the high 
H end. At this point we believe this should be regarded as 
an empirical question that can only be answered by future 
high-resolution studies with instruments such as ALMA. 

Another potential source of systematic uncertainty lies 
in the fact that there is only one Millennium Simulation, and 
it was not carried out for the currently best-fitting cosmo- 
logical model. In particular the value of as used (the normal- 
ization at a scale of 8 h~ Mpc) was 0.9, wh ile the WMAP5 
value is more like 0.8 (|Dunklev et alj |2008). A smaller nor- 
malization can have a significa nt effect in reducing the in- 
ciden ce of strong lensing (e.g., iLi et al . 2007; Fedcli et al.l 
2008). However, it is difficult to assess precisely what affect 
this might have on our estimates, for several reasons. Firstly, 
the full simulations are not available for other models. Sec- 
ondly, there are other differences in the current best-fitting 
model, e.g. Qm and the primordial fluctuation slope, n (a red 
tilt means there are lower fluctuations at galaxy and group 
scales than in the Millennium simulation). And finally, there 
is still some uncertainty about the precise value of as, with 
the statistics of gian t arcs suggesting that higher valu es may 
be a better fit (e.g. jLi et al.ll2007l ; iFedeli et al.ll200Sl ). 
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5.2 Source counts 

As already discussed, we have taken the shallowest power- 
law consistent wit h the SHADES data as found in 
ICoppin et al.l l|2006l ) and used that as an upper estimate on 
the source counts at the bright end. At high flux densities, 
this upper estimate and our best guess Schechter function 
give significantly different results, with the 'maximal' counts 
providing the lower estimate for lensing effects. 

It is possible that the observed counts already include 
some magnification bias. This being the case, our estimation 
of the expected number of bright sources will be too high, 
since these already amplified objects will be lensed further 
by our procedure in Section [4] To attempt to compensate 
for this, we looked at how the slope of the differential num- 
ber counts changes when lensing is applied only at the flux 
densities where the data exist. Using the Schechter function, 
we found a negligible change, from which we conclude that 
lensing effects in this region are too small to change our re- 
sults. However, for the broken power-law, the slope after the 
break (/3 in equation [TJl is much more sensitive to changes 
introduced by lensing. Above the knee, ICoppin et all (|2006h 
found [3 — 5.1 ± 0.9. Fitting equationQto the lensed source 
counts for Sobs < 25, where the counts still follow a broken 
power-law and SHADES data exist, we find f3 = 4.9. This is 
still well within the uncertainty in the original broken power- 
law. In p articular, the limit of 7V(>22mJy) < I7deg -2 
found by ICoppin et al.l (|200rJ ) is maintained by the maxi- 
mal broken power-law even after lensing. 

We have been making the simplifying assumption that 
the shape of the submm counts is independent of redshift, 
even though the amount of lensing might follow a particular 
source redshift distribution. The main difficulty in carrying 
out estimates within more complicated models is that there 
are only relatively weak constraints on how the shape of the 
counts might vary with redshift. The best available evidence 
comes from a study of lumino sity evolution among SMG s 
in the GOODS-North field by I Wall. Pope fc Scottj [|2008h . 
which builds on earlier suggestions that the brightest SMGs 
tend to be at higher redshift. This appears to run counter to 
the idea of hierarchical clustering, except of course it is most 
likely complications with the baryon s within ga l axy fo rma- 
tion that lead to this 'down-sizing'. IWall et all <|2008l ) find 
evidence for 2 populations separated in luminosity, each hav- 
ing a different luminosity function slope and redshift depen- 
dence. To the extent that the negative if-correction makes 
up for distance dimming for SMGs, the submm flux den- 
sity is a proxy for luminosity, with the break occurring near 
5 mjy: fainter SMGs have a shallower slope and evolve more 
gently; brighter SMGs have a steeper slope and evolve more 
strongly with redshift. Using thi s approximation w e find that 
the phenomenological models of lWall et al.l (|2008h do indeed 
lead to SMG counts which vary with redshift. However, the 
results are (within the fairly wide error bars) consistent with 
our 'maximal' model. 

At the moment the available SMG data involve modest 
sample sizes, limited dynamic range and incomplete redshift 
information. As more is learned about how the counts evolve 
with redshift it will be possible to carry out more accurate 
calculations of the lensing effects. 



Q_ 




10 20 30 40 50 60 70 80 90 100 
Observed flux density S obs [mjy] 



Figure 10. Dashed lines are the fraction of dN /dS b s due to 
magnifications greater than 2 using the model from Appendix [A] 
in redshift bins of Az = 0.5 for 0.5 < z < 5.0. Solid lines are from 
the same redshift intervals using the Schechter function. The thick 
lines are the average over all redshifts. 



6 COMPARISON WITH AN EVOLUTIONARY 
MODEL 

A short-coming of our analysis is the assumption of a fixed 
shape for the source counts at all redshifts. Naive expecta- 
tions for galaxy formation are that the most extreme SMGs 
would be rarer at high redshift. However, the idea of 'down- 
sizing' has been shown to apply to SMGs (Wall, Pope & 
Scott 2008), making the 'hyper-LIRGs' more prevalent at 
earlier times. Hence, to properly assess these effects requires 
a more detailed model which includes the evolution of SMGs. 
We choose a particular set of simulations based on the lo- 
cal far-infrared colour-luminosity distribution described in 
Chapin, Hughes & Aretxaga (2009) including evolution that 
reproduces the observed 850 u.m source counts and redshift 
distribution (see the complete description in Appendix |A")) . 
However, since the existing empirical information is so mea- 
ger, the model's behaviour at high redshift can only be con- 
sidered as indicative. Moreover, since the model is based on 
a particular simulation, then it is limited for rare objects 
at the bright end. Nevertheless, its use should enable us to 
assess the effect on our estimates of allowing the number 
count shape to vary as a function of redshift in a way which 
is consistent with current data. 

The simulation we utilise contains over 8 million sources 
between z — and 5. We divided these sources into bins 
of width Az — 0.5 and determined the differential source 
counts dN/dS for each bin, to which we could apply the 
lensing probabilities in the same way as before. The prob- 
ability of a lensing event as a function of Sobs is shown in 
Fig. 1101 and compared to the Schechter function using the 
same redshift bins. 

It can be noted that the evolutionary model produces 
considerably more variation in lensing probability compared 
to the analytical Schechter function. Nonetheless, the prob- 
ability is still greater than about 0.3 for any object with 
an observed flux greater than about 20mJy, which is com- 
parable to the results we found for a single counts shape. 
For dimmer objects, where we may need to be conservative 
about what sources deserve closer attention, the chances of a 
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lensing event are actually slightly better with the evolution- 
ary model. So the conclusion is that although our analytical 
approximation misses some of the structure in the redshift 
distribution of SMGs, the end results are quite similar. Still, 
this should be considered as yet another potential source of 
uncertainty in the lensing predictions, and further motiva- 
tion for gathering data which can distinguish among details 
of the models. 



7 DISCUSSION AND CONCLUSIONS 

We can use our study to make approximate predictions for 
upcoming submm surveys. There are several relevant in- 
struments, but we focus on the s urveys which a r e plan ned 
with the SCUBA-2 instrument ^Holland et all l200rj ) on 
the James Clerk Maxwell Telescope. There are 2 relevant 
850 /J,m programmes: the SCUBA-2 Cosmology Legacy Sur- 
yey, S2CLS; and the SCUBA-2 'All Sky' Survey, SASSy 
ijThompson et alJl2007t >. 

We have also carried out estimates for surveys at shorter 
wavelengths, of the sort which might be performed with the 
450 /um array of SCUBA-2, or the SPIRE instrument on the 
Herschel satellite (operating at 250, 350 and 500 pirn). While 
it is clear that there may be many examples of strong lenses 
in such wide surveys, the fraction of bright sources which 
are lensed is significantly lower than at longer wavelengths. 
If one wants to find such 'monsters', either as probes of 
line-of-sight structure or for their own intrinsic value, then 
one should turn to ground-based surveys in the 850 /jm or 
~ 1 mm windows. 

The S2CLS plans to map approximately 20deg 2 to an 
RMS of 0.7 mjy. From our number counts model we esti- 
mate that there will be about 96, 44 and 27 sources de- 
tected with S > 20, 25 and 30 mjy, respectively. The total 
number of sources above these flux limits which have n > 2 
will be about 15 (with most of them at the high flux end) 
using our best estimate from the Schechter function num- 
ber counts. The numbers change only slightly over this flux 
density range, due to the probability of lensing increasing 
faster than the decrease in source counts. Above 35 mjy the 
number of expected lensed sources declines steadily. The 
lensed fractions are considerably smaller using our 'maxi- 
mal' counts model, as can be seen by referring to Fig. [6] 

There will of course always be bright sources which have 
negligible lensing - hence one would like to know how bright 
to go before the probability of strong lensing is significant. 
This can be determined using Fig. [6] If one is prepared to 
accept a 1 in 3 chance of a source being strongly lensed, then 
one should select sources observed with S > 25 mjy. If one 
would like the chance to be 1 in 2, then that flux rises to 
about 30mJy. 

What would one do with such sources in practice? Since 
strong lensing is likely to come from either a galaxy clus- 
ter, or a massive galaxy, then the existence of strong lens- 
ing implies strongly clustered structure along that line of 
sight. Hence for each sufficiently bright source one would 
use follow-up observations at other wavelengths to try to 
establish whether strong lensing was likely, and then to see 
if one could find the structure which was responsible for the 
lensing. Multiple images or distorted morphologies of optical 
counterparts would be ways of determining that lensing was 



taking place - these would naturally show up as part of the 
procedure for trying to determine counterparts in deep data 
at other wavelengths. For SMGs where lensing was strongly 
suspected, one would target the area to search for the pres- 
ence of structure along the line of sight, either with X-ray, 
Sunyaev-Zel'dovich or 'red cluster sequence' observations. 
The advantage of this approach is that the high redshift of 
the SMC sources means that it should be feasible to find 
cluster (or proto-cluster) lenses at higher redshifts than are 
easy to achieve with most other techniques. Of course, the 
selection effects for clusters found in this way may be compli- 
cated to quantify. Nevertheless, building up samples of z > 1 
clusters is sufficiently important for understanding structure 
formation (as well as constraining dark energy, etc.), that it 
is worth using every available method. 

SASSy is designed to make a shallow 850 /im map over 
approximately 4,000 deg 2 , or one tenth of the sky, with the 
possibility of extension to a larger area later. The RMS of 
the maps is planned to be around 30 mjy, so that a robust 5a 
catalogue will have a limit around 150 mjy. It may also be 
possible to reduce this to nearer to 100 mjy using targetted 
repeat observations for peaks in the maps. 

The procedure for identifying strongly lensed sources in 
SASSy may be a little different than for S2CLS, and be- 
cause of the brighter flux densities, the level of uncertainty 
in predictions of the number of lensed sources will be con- 
siderably higher. These predictions depend strongly on the 
counts model, the normalization of the high amplification 
tail of the lensing PDF, as well as the amplification cut-off 
imposed by the finite source size for SMGs. Since there are 
still huge uncertainties in all of these factors, the expecta- 
tions for SASSy cover a wide range of possibilities. 

Using the optimistic limit of 100 mjy for SASSy, the un- 
lensed model counts give approximately 1200 sources for the 
SASSy catalogue. Many of these will be in the 'Euclidean 
counts' regime, and hence should be relatively easy to elimi- 
nate as lensing candidates. These will often be already well- 
known galaxies, with others typically being in the IRAS, 
Akari or radio catalogues. Colours (e.g., 850 /im to radio) 
can be used to distinguish objects which are likely to be 
at higher redshift and hence have a higher likelihood of be- 
ing lensed; such methods will also be necessary to eliminate 
Galactic clouds. Our best lensing estimate yields a total of 
1000 sources in SASSy with lensing amplification n>2, with 
most of them being much more strongly lensed than this 
limit. We expect that these extremely lensed sources will 
be fairly easy to distinguish from relatively nearby, intrinsi- 
cally bright galaxies, and hence the chances of there being 
structure along the line of sight to such candidates will be 
very high. Follow-up of these SMGs at other wavelengths 
will also be easy, since they should be at least an order of 
magnitude brighter than the typical SCUBA sources which 
have been followed up in the past. 
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APPENDIX A: EVOLUTIONARY MODEL 

The evolutionary model used to determine 850 [im counts as 
a function of redshift combines the local measurement of the 
joint far-IR colour-luminosity distribution, $(L, C) (L is the 
42.5-122.5 um far-IR luminosity, and C is th e logarithm of 
the 60 /100 um flux density ratio) presented in I Chapin et al.l 
(2009), with a two-popul ation evolutionary scenario similar 
in style to the model of lLagache et al.l l|2003l . |2004f) , con- 
sisting of lower-luminosity 'normal' galaxies, and higher- 
luminosity star-bursting galaxies. The constraints for thi s 
model are the sou rce counts at 850 um jCoppin et alj feooe), 
and 170/160 um (|Dole et all 1200 ll . |2004|), and the redshif t 
distributions of 850 um sources from Chapman et alj (2005). 
In addition we verify that the integrated contribution to the 
IR background in the ra nge 850-160 um do es not exceed 
that measured by C OBE dFixsen et al.|[i99ST ). 

First, similar to lYun et al.l (|200 ll ). we approximate the 
local far-IR luminosity function by the superposition of two 
Schechter functions, (j>(L) = p*(L/L*) a exp(— L/L t ). For 
the lower-luminosity objects p, = 2.3 x 10~ 3 Mpc -3 dex -1 
and L* = 2.3 x 1O 1O L0, while for the higher-luminosity 



objects p, = 2.3 x 10" s Mpc" a dex" 1 and L, = 2.2 x 
10 11 L©. A value of a — —0.5 is adopted for both popula- 
tions. T he two compo n ents a re compared with the original 
form of IChapin et ail (|2009T ) in Figure [M] The full local 
colour-luminosity distribution for each of these populations, 
C) and &2(L, C) are then constructed by multiply- 
ing the respective Schechter functions by the conditional 
co lour distribut i on p(C \L), as defined in Equations 8-10 
in I Chapin et all i|2009h and using coefficients for the 'dual 
power-law luminosity function' fit. 

The 'normal' galaxy population is evolved in redshift 
using a combination of density and luminosity evolution, 
$i(L,C, z) = g(z)<&2[L/h(z), C]. For the density evolution 
we use the form g(z) = (1 + 2) 3//2 sech 2 [&ln(l + z) — c] cosh 2 c 




Figur e Al. The local 42 .5-122.5 pm far-IR luminosity function 
from IChapin et al . (2009) expressed as the superposition of a 
low-luminosity 'normal' galaxy population, </>i(L), and a high- 
luminosity starburst population, ij>2(L), each parameterized with 
Schechter functions. 
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Figure A2. Observe d integral sour ce counts at 850 pm (triangles 
with error bars.lCoppin et all200 6h and a t 170/160 (am (triangles 
without error bars. lDole et alj20 01. 2004). The model predictions 
are given by the solid lines. For the purposes of this comparison 
the 170 um data are also assumed to be at 160 (jm (which is a 
reasonable approximation given the proximity in wavelength). 



from iLewis et all l|2005h . since the 2-parameter family of 
curves smoothly describes a wide range of scenarios involv- 
ing a steep rise with a gradual fall. We adopt coefficients 
b — 4.5 and c = 1.6. For the luminosity evolution h(z) grows 
as (1 + z)' 2 ' 5 up to z — 2.0, remains constant to z = 4.0, and 
there are no more galaxies at higher redshifts. The starburst 
galaxy population is defined using pure luminosity evolution, 
$2(£, C, z) — <&2[L/g(z), C], where g(z) is the same param- 
eterization used for density evolution of the low-luminosity 
population, with b = 1.5 and c = 1.3. 

Within this model the colour-luminosity correlation un- 
dergoes the same luminosity evolution as t he luminosity 
functi ons for each population. It was shown in lChapin et al.l 
(2009) that such evolution can explain the observed temper- 
atures for a sample of 850 um-selected galaxies at redshifts 
z> 1. 
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Figure A3. The observed peak-normalized redshift distribu- 
tion of 850 um sources corr ected for selection effects (histogram) 
from lChapman et al.l d2005T) , compared with the model prediction 
(solid line) for sources with 5850 > 5 mjy. This flux limit approx- 
imates that of the SCUBA surveys from which the sources were 
selected for the redshift survey. 



SED templates are required to calculate flux densi- 
ties for particular galaxies with luminosities L and far-IR 
colours C drawn from the evolving distribution. Since we 
only consider data in the oberved wavelength range 850- 
160 urn we are free to us a simple modified blackbody SED, 
S v ot u^B^lT), as there is no significant influence on the 
shape from hotter dust components, or complicated features 
related to AGN activity or p olycyclic arom a tic hy drocarbons 
in the mid-infrared. As in IChapin et al.l (|2009) we adopt 
f3 = 1.5 such that C maps directly to an effective dust tem- 
perature T. 

The model predictions for 850 and 160 um source counts 
are shown in Figure lA"2l along with existing measurements, 
demonstrating the agreement across this broad range in 
wavelengths. In Figure IA3I we compare the model redshift 
distribution for 850 urn-selected sou rces brighter than 5 mj y 
compared with the redshift survey of lChapman et all l|2005h , 
also exhibiting good agreement. 

With this model set up to be consistent with current 
knowledge of the local luminosity function, submm number 
counts and redshift distributions, we can use a realization of 
it to generate a simulated source catalogue for studying the 
effects of redshift evolution on the lensing statistics. 

This paper has been typeset from a TgX/ WFpfi file prepared 
by the author. 



